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Artificial cavity photon resonators with ultrastrong light-matter interactions 
are attracting interest both in semiconductor and superconducting systems, 
due to the possibility of manipulating the cavity quantum electrodynamic ground 
state with controllable physical properties. We report here experiments show- 
ing ultrastrong light-matter coupling in a terahertz metamaterial where the 
cyclotron transition of a high mobility two-dimensional electron gas is cou- 
pled to the photonic modes of an array of electronic split-ring resonators. We 
observe a normalized coupling ratio ^ = 0.58 between the vacuum Rabi fre- 
quency Q and the cyclotron frequency uu c . Our system appears to be scalable 
in frequency and could be brought to the microwave spectral range with the 
potential of strongly controlling the magnetotransport properties of a high- 
mobility 2DEG. 
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Enhancement and tunability of light-matter interaction is crucial for fundamental studies 
of cavity quantum electrodynamics (QED) and for applications in classical and quantum de- 
vices El IH IH HI . The coupling between one cavity photon and one elementary electronic 
excitation is quantified by the vacuum Rabi frequency f2. The non-perturbative strong light- 
matter coupling regime is achieved when is larger than the loss rates of the photons and elec- 
tronic excitations. Recently, growing interest has been generated by the ultrastrong coupling 
regimeflU [H HI [TOl EEU CGI, which is obtained when the vacuum Rabi frequency becomes 
an appreciable fraction of the unperturbed frequency of the system lu. In such a regime, it is 
possible to modify the ground and excited state properties obtaining non-adiabatic cavity QED 
effects 0. Experimental progress has been achieved in two different solid-state systems : (i) 
microcavities embedding doped quantum wells lfP3l IT4l [T5l [T6l [T71 . where the active electronic 
transition is between quantized subbands in the well; (ii) superconducting quantum circuits in 
transmission line resonators lfT8l fT9l , where the photon field is coupled to artificial two-level 
atoms obtained with Josephson junctions. 

We present experimental results on a new system, namely a high-mobility two-dimensional 
electron gas (2DEG) coupled to terahertz (THz) metamaterial resonators. The photon mode 
is coupled to the magnetic cyclotron transition of the 2DEG, obtained by applying a magnetic 
field perpendicular to the plane of the quantum wells (Fig. 1(A)). The cyclotron frequency is 
expressed by co c = ^ where B is the applied magnetic field, e is the elementary charge and 
m* represents the electron effective mass. This highly controllable system is ideal for the study 
of strong coupling because the material excitation can be continuously tuned by changing the 
value of the applied magnetic field. The key physical aspect to highlight is the dependence of 
the optical dipole moment d for a cyclotron transition on the cyclotron orbit length. The dipole d 
scales as d ~ eloy/u [20], where l = ^Ti/eB is the magnetic length and v = P2Deg2vt/q is the 
filling factor of the 2DEG , being P2DEG the electron areal density. This proportionality of the 
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dipole with respect to Zo allows to have gigantic dipole moments as soon as the cyclotron tran- 
sition can be resolved. According to theoretical calculations valid for integer filling factors and 
for an optimized resonator geometry, the coupling ratio is expected to scale as ^- ~ ^Janq^u 
where a is the fine structure constant and ?t,qw is the number of 2DEGs [20J. For high filling 
factors, this coupling ratio is predicted to assume values even larger than unity (corresponding 
to transitions in the microwave range). We use high-mobility 2DEG based on GaAs material 
system [|2T| and we realize our experiments in the THz region of the electromagnetic spectrum. 
These frequencies, for our material system, correspond to magnetic fields of the order of a few 
Tesla and optical experiments are conducted using broadband THz pulses generated with ultra- 
fast lasers (221. A THz-TDS system (bandwidth 0.1 - 3 THz) [23] is coupled to a split-coil 
superconducting magnet to probe sample transmission [|24l . 

Our THz metamaterial integrates the 2DEG with a metasurface of electronic split-ring res- 
onators (see Fig. 1(A)). These resonators [|26l 1271 |28l exhibit electric field enhancement over 
strongly subwavelength volumes [29], making them ideal candidates to reach extreme couplings 
in the Mid-Ir and THz range where long wavelength radiation has to interact with quantum well 
systems typically extending over length of some micrometers [30]. Moreover, the enhanced in- 
plane (x-y) electric field couples efficiently to the cyclotron transition when the magnetic field 
is applied perpendicularly to the plane of the layers and parallel to the wavevector of the inci- 
dent THz pulse (see Fig. 1(A)). Resonators were deposited on top of the 2DEG by conventional 
photolitography, metallization with Ti/Au (5/250 nm) and lift-off technique. 

At zero magnetic field we observe two resonances m l5 m 2 whose origin is qualitatively 
different: the lowest frequency mode {fx m 0.9 THz ) is attributed to the LC resonance, where 
counterpropagating currents circulate in the inductive part and the electric field is enhanced 
mainly in the capacitor gap [27]. The second mode (f 2 ~ 2.3 THz) is attributed to the "cut 
wire" behavior where a A/2 kind of resonance is excited along the sides of the metaparticle 
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H28B . These values correspond well to simulations with 3D FE modeling (see SI in H24I0 . The 
two modes of the split-ring resonator also have different transverse wavevectors, as it is evident 
looking at the different field distributions. The presence of conductive layers underneath alters 
the frequency and the quality factor of these resonances. 

We observe a value of QrriL — ^-3 for an insulating substrate and Q^^ 1 ~ 3.1 when the 
resonator is deposited on top of the single 2DEG sample. The values for the second resonance 
results less affected yielding Q^jthz — 5-3 and Q^stuz — 5-3 (a more detailed analysis can 
be found in SI of 11241 ). It is important to highlight that, in contrast to atomic systems, we can 
realize strong light-matter coupling physics with resonators displaying extremely low quality 
factors. The giant value of the coupling constant Q ~ v MqwP2Deg typical of intersubband 
systems [5| together with the high electric field enhancement of sub- wavelength metallic res- 
onators (V cav ~ 8 x 10" 17 m 3 in our case ) allow the observation of cavity polaritons in a system 
where both components are in principle highly dissipative. 

In the data reported in Fig. 2(A) we observe the evolution of the sample transmission \t\ = 
| e^oeg^o) I as a f unc ti° n °f me applied magnetic field (normalized to the electric field £^2Deg(0) 
of the reference 2DEG wafer at B = 0T) . One 2DEG (nQw = 1, doping density P2DEG = 
3.2 x 10 11 cm -2 ) is used as an active medium and placed 100 nm below the surface: its cy- 
clotron resonance can couple to the resonator modes. As the magnetic field is swept a profound 
modification of the sample transmission is observed. The possibility to tune in a continuous 
way the material excitation allows to follow the evolution of polaritonic states as the system is 
driven from the uncoupled regime to the strongly coupled one [13111 . We observe two successive 
anticrossings when the cyclotron energy matches the first and the second resonator modes. 

In Fig. 2(B) we extract the positions of the minima of sample transmission and plot the 
dispersion curves for the polariton eigenvalues as a function of the magnetic field. The curves 
are calculated using a full quantum mechanical treatment of the system, obtained generaliz- 
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ing the theory described in Ref. ll20ll to the case of a zero-dimensional resonator exhibiting two 
modes with different transverse wavevectors |[24]|. Following a bosonization procedure, we have 
derived the different contributions to the total Hamiltonian and have diagonalized it using the 
Hopfield-Bogoliubov method. It is worth mentioning that the ideal resonator we have consid- 
ered in the analytical treatment is different from the real split-ring one. However, we emphasize 
that this would introduce only a form factor in the matrix element calculation. 

In order to fit the experimental data, we need to know the resonator modes frequencies (/{ 
and f 2 ) as well as the strength of their couplings (f^ and VL 2 ). For each cavity mode we 
assumed the asymptotic value of the corresponding lower polariton branch to coincide with 
the frequency of the unloaded resonator (f[ = 0.83 THz and f 2 = 2.26 THz) flMEU • The 
coupling strength - for the two modes can not be directly measured, we thus applied a best 
fit procedure following the least square method, analogously to what done in Ref. Ifl4ll (more 
details in S3 of [24]). The minimal error O is obtained for Oi = 0.17 and Oa = 0.075 
(where to-y = 2nf[ and w 2 = Zirfy)- As expected the coupling strength scales monotonously 
with v : for the measured density P2DEG' we nave V {E) = v(2 T) ~ 6.5 for the first mode and 
v(B) = z/(5.5T) ~ 2.4 for the second one. 

To increase the coupling strength, we kept the resonator geometry and hence the frequency 
constant, and increased the effective number of carriers in the system. A new sample was 
prepared with r2 QW = 4 wells and a doping density per well /4deg = 4-45 x 10 11 cm~ 2 (see 
scheme in Fig. 1(B) and materials and methods in ll24l ). Sample transmission as a function of 
the applied magnetic field (Fig. 3(A)) shows that the system is driven deeply into the ultrastrong 
coupling regime. The polaritonic line widths display narrowing as the low quality cavity mode 
is mixed with the cyclotron resonance (Fig. 3(B)). Following the fitting procedure previously 
described we observe a value of — = 0.36 for the first resonance. Indeed, the effects of the anti- 
resonant terms of the light-matter Hamiltonian start becoming relevant when the dimensionless 
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ratio - is of the order of 0.1[14]. Due to the increased doping, the filling factor in the region 



2.11 and theoretically J (° EG = 2.35. The small discrepancy can be attributed to the different 



coupling of the quantum wells which do not experience the same electric field of the resonator: 
the generalized expression of the coupling ratio is calculated in the case when all the wells are 
coupled in the same way to the resonator's field Il20ll . 

By further scaling the resonator frequency down to / ~ 500 GHz employing a slightly 
modified geometry (see inset Fig. 3(D), section S2 and Fig. S5 of ||24) ) and employing the 
sample with r\q W = 4 quantum wells we could probe the regime where the polariton splitting 
at the anticrossing 2hVt is larger than the bare cavity photon energy. In Fig. 3(D) we report the 
positions of the minima of the sample transmission for the case of f=500 GHz resonator together 
with the fitted dispersion curves: we measure a normalized ratio ^ = 0.58 for a filling factor 
of v{1.2T) ~ 15.2 which corresponds to 2htt ~ 1.2w c (see also Fig.S6 in [EH ). 

The generalization of the theory developed in Ref. Il20l accounts for the depolarization 
shift in presence of a magnetic field (magnetoplasmon) originating from the long-wavelength 
part of the Coulomb interaction. We found that the renormalization of the cyclotron transi- 
tion frequency is too small to allow the experimental resolution of the magnetoplasmon branch 
(see Figs 2(A) and 3(A)) in our experimental parameter regime including the small wavevector 
condition ql <C 1 which is always satisfied since dealing with optical wavevectors. 

We have observed ultrastrong light matter coupling in a composite THz metamaterial mea- 
suring a normalized coupling ratio ^ = 0.58. The impact of our results has to be considered 
also in the perspective of a change in the DC transport properties of the 2DEG, in analogy with 
what already observed by direct irradiation at lower frequencies [l32l . 



of the anti crossing is v ~ 9. As expected, the coupling ratio scales with a /P2Deg^qw and for 



the two samples at the same resonant frequency we obtain experimentally 
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These results should lead to the scaling of the frequency to lower values and to an increase 
of effective density to further enhance the coupling strength. 
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Figure 1: (A): Schematic of the composite metamaterial used in our experiment together with 
the experimental arrangement showing the polarizing static magnetic field, the wavevector and 
the polarization of the incident broadband THz pulse, the Landau level arrangement and the 
semiclassical representation of a cyclotron orbit of magnetic length l . A metasurface com- 
posed by LC-metaparticles with a design similar to Ref.[25] is deposited on top of the semi- 
conductor. SEM picture of one metaparticle: the split gap of the capacitance elements is 
2.6 fim. (B): The bandstructure of the multi 2DEG system is schematized together with the 
quantum well position (not to scale) and (C) x-y spatial distribution of the in-plane electric field 
(E P iane=\J\E x \ 2 + \E y \ 2 ) calculated with a finite element commercial software is shown for the 
two observed modes mi and m 2 . (z=100 nm below the semiconductor surface surface), (d): 
Intensity for E p i ane in the yz plane for the low frequency mode mi (cut along the white dashed 
line in (C) mi ). 
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Scalari, fig2, x43p4 




Magnetic Field [T] Magnetic Field [T] 

Figure 2: ((A): Transmission \t\ of the sample («qw = 1) as a function of the applied magnetic 
field B. The reference is a plain 2DEG sample without resonators on top and the measurement 
is performed at T=2.2 K. (B): Best fit with the extracted transmitted minima positions for the 
two different transverse modes of the electronic split ring resonator; the fitting parameter is — . 
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Figure 3: (A): Transmission \t\ of the sample (nqw = 4) as a function of the applied magnetic 
field. The reference is a plain 2DEG sample without resonators on top and the measurement is 
performed at T=10 K. The black dotted line highlights the cyclotron signal coming from the un- 
coupled material which is left in-between the resonators. (B): Sections in the two anticrossing 
regions for the sample transmission. (C): Best fit with the extracted transmitted minima posi- 
tions for the two orthogonal modes of the split ring electronic resonator; the fitting parameter 
is — . (D): Best fit with the extracted transmitted minima positions for the f=500 GHz resonator 
and n QW = 4 measured at T=10 K ; the fitting parameter is Inset: scheme of the 500 GHz 
resonator. 
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